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are given in genus zero by Legendre polynomials whose argument is determined by 
an analog of the string equation. 
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1 Introduction 



Recently the external field problem for matrix models has achieved a remarkable interest 
due to the work by Kontsevich who had represented the partition function of 2D 
topological gravity as that of a hermitian one-matrix model in an external field. Along this 
line, a direct proof of the Witten's conjecture about an equivalence of 2D topological 
and quantum gravities has been obtained 0, H, ||. 

On the other hand, the Kontsevich matrix model can be explicitly solved |^, |^ 
order by order of genus expansion making use of the standard methods which were first 
introduced for an analogous unitary matrix problem M. This solution has been applied 



recently by Itzykson and Zuber to calculate explicitly the intersection indices on the 
moduli space. 

The original Kontsevich model is reduced by a linear shift of an integration variable 
to a generic form 



^[A]=|DXe^*'-(^^+^™ (1.1) 



with a cubic potential V{X) and the integral going over N x N hermitian matrices X. For 
an arbitrary V{X), the partition function ( |1.1| ) is associated with the so-called generalized 
Kontsevich model. As is advocated by Kharchev et al. this model interpolates 

between arbitrary [p, q) double-scaling limits of the standard multi-matrix model. It is 
crucial for applications that the exponent in (|1.1|) depends on the external field A linearly. 

One might think, however, about another way of incorporating the external field by 
modifying the quadratic term: 

Z[A] = j /)Xe^*'-(-^™+^(^)) . (1.2) 

A model of this kind was first introduced by Das et al. [O. In the language of discretized 



random surfaces, it corresponds to an external field which is connected to curvature of 
the world sheet. 



One may ask a question: which form of the potential V{X) in ( |1.2| ) permits such 
changing the variables that results in an integral of the type ( |1 . 1| ) ? The only nontrivial 
choice is the nonpolynomial V{X) = a(log(l + X) — X) which is associated with the 



Penner model |T3[. As is pointed out by Distler and Vafa |T^, the double-scaling limit of 
the Penner model corresponds to c = 1 CFT or to the d= \ string at the self dual radius 
R = \ which is first solved by Gross and Klebanov |T^. Further studies of generalizations 



of the Penner model at multi-critical points have been done recently in |]T6|, |17|. 

We solve the model ( [L.2D in genus zero explicitly using the method of Brezin and Gross 
0. We show that the second derivative of log ^ in a is a pure logarithm which is a feature 
of c = 1 theories. We study various multi-critical regimes and relate the traces of inverse 
powers of A (the Kontsevich-Miwa variables) to sources for conformal operators. We 
calculate the correlation functions of the conformal operators with the puncture operator 
and show that it is given by Legendre polynomials whose argument is determined by an 
analog of the string equation which is in our case a set of two equations for two variables. 
We calculate the proper anomalous dimensions which turn out to be different at different 
multi-critical points. 
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2 Formulation of the model 



We consider the hermitian one-matrix model in an external field: 

Z[A]= DXe \' L 81 J J/^ A = diag(Ai,... ,A^). (2.1) 

Our strategy to deal with the model (^TTj) is to reduce it to the Kontsevich type integral 
(pni). This can be done by a substitution A^/^XA^/^ — > X. We obtain 



Z[A] = e-Tf-A'+"^'(detA)-^(°+i) X 



(2.2) 



Now our goal is to study the integral 

Z[r]\ = J Me-'^*''{5^'-''^}(detX)'^^. (2.3) 

There are different techniques to proceed with this integral. In the next sections we shall 
treat it using the method of Schwinger-Dyson equations written in terms of eigenvalus, 
which has been applied to hermitian one- matrix models in p, ^ , and discuss Virasoro 
constraints for this model. 



We conclude this section with the Itzykson-Zuber-Mehta technique for integration 
over angular variables in multi-matrix models. We assume that it is well known and do 
not intend to engage into any detailed description. In terms of eigenvalues of the matrices 
X and T], (|2.3|) can be expressed as follows: 

m = Z rfx.^M exp |-Ar;£ - ry,x.) | , (2.4) 

where A [a:] = ni>, ( 

the Van der Monde determinant. In what follows we assume 
that aX is positive integer. Using the relation 

r dx 6-"'/'+''" = r" Hn{i7])e''''^ (2.5) 

J —oo 

we may integrate out all Xj's in the Eq. (|2.4| ): 



N ^2 

l<j<N ^ i=l 



Z[r]] = 2iv°^V2+JV(A^+i)/2^-i[^] ^^^^^ ||^-"^-('=-i)if,;v+fc-i(^v^r/i)||exp(^iV5^| 



(2.6) 

Here Hn{ii]) are the Hermite polynomials of an imaginary argument. As is shown by 
Kharchev et al. [11|, the determinant formula ( p.6|) implies Z[r]] to be a r-function. 



It is clear from Eq. ( p.6|) that Z[ri] is a symmetric polynomial of degree aN"^ in variables 
rji. But it is still unclear how to deal further with the determinant. More effective way to 
obtain the answer is to use the Schwinger-Dyson equations generating Virasoro constraints 
and the "master equation" in this model. 
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3 Virasoro constraints 



The Schwinger-Dyson equation for our model ( |2.1|) follows from an invariance of DX 
under (X-independent) variations of X: 



f 5 -N tr(^-r)X-a\ogX] / 6 f X^ , ^'^ 



(3.1) 



Taking into account that < XijF >= < F > and applying one additional 
external derivative over 77 to eliminate the nonlocal contribution arising from the variation 
of alogX, we have 



o2 / \ 



Z[r,] = 0. (3.2) 



This equation is discussed in detail below. 

One could make alternatively the X-dependent variation 5X = enX""*"^ which would 
immediately lead for n > to a set of Virasoro constraints of the type advocated by 
Semenoff and one of the authors for the model ( p. . 1| ) with arbitrary polynomial V{X): 

CnZ[ri] = for n>0 (3.3) 

and 

= E I -(V.)"+' + N{V,r+% + aN{V,r + ^ E E(VO'(V,)"-M (3.4) 

where 

V. = ^ + E^— • (3-5) 

9Vt r]i - r]j 

However, the n = —1 equation which corresponds in the polynomial case to £_i operator 
is now nonlocal. This seems to be a reflection of the fact that we are dealing with a c 7^ 
theory. 

Let us now turn to Eq. (|3.2| ) which resembles in many details the equation 
[S"^ /Sr]"^ + 7]) Z[ri] = generating |^ the continuum Virasoro constraints in the stan- 
dard Kontsevich model. For our case we do not know in advance proper expressions for 
continuum time-variables in terms of the eigenvalues of A. We show below by explicitly 
solving the model ( p.l|) in genus zero that those are given by the Kontsevich-Miwa trans- 
formation while the problem of constructing continuum Virasoro (or W-type) constraints 
for our model is beyond the scope of present publication. 

However, it is worth mentioning that for the model (|2.1| ) there exists some Virasoro 
algebra in terms of "reverse times" Qk = ^ trt]^ while the standard Kontsevich-Miwa 
times are defined via expansion in negative powers of the matrix 77. 

Doing the contraction of Eq.( p.2| ) with {ri"')ki we obtain the set of constraints: 

00 

J2<ls+n~2LsZ[7]], (3.6) 
s=2 
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where 

r\'2, OO ^1 ^1 

a+b=s (^(laOQb Clqs+b OQs Oqs-2 

One may prove that for finite these conditions are sufficient to fix uniquely the form 
of the solution which is the polynomial of the maximum order aiV^ in qi multiplied by 
a trivial Gaussian factor. The continuous limit is rather nontrivial and will be discussed 
elsewhere. 

We shall consider in what follows not the Eq.( |3.2|) but the equation for the ratio of the 
determinants ( p.6| ). Z[ri] differs from it by a factor e~ ^^^^ . Pulling this factor through 
Eq.( |3.2|) , we have 

-i^] -N(v-f] +aN\, 



det 



HaN+k-iiWNr] 
""AM 



0. 



There also exists an interesting representation of these constraints not in Kontsevich- 
Miwa but rather Macdonald's like variables [|I9|] : 

9k = Vh---Vik, (3-9) 



where the sum runs over all products of k different r/j's. Then the symmetric function in 
Eq.( |3.8| ) is the polynomial of degree aN (not aN"^) in the variables Qk, the higher term 
being gf/^ . Moreover, the coefficient standing by term QN-ki ■ ■ ■ dN-k^M depends only on 
the set {ki, . . . , kaN} but not only itself. In this case the exphcit solutions for few 
lowest values of aN can be found. 



4 The "master equation" 



In this section we rewrite Eq.( ^.2[ ) (or ( ^.8D ) in terms of eigenvalues of the matrix i]. We 
obtain an equation (the "master equation" ) which is similar to that of Refs. [0, Q for 
the Kontsevich model. 

We want to find a suitable form of Eq. (|3.8| ) in terms of eigenvalues of the matrix rj: 

r] = UHU+, if = diag (r/i,...,r7^), (4.1) 

where ?7 is a unitary matrix. We contract Eq.( p.8| ) with the special matrix F^^^ of the 
form: 

FO) = f//,f/+, /, = diag (0, . . . , 0, 1(,), 0, . . . , 0). (4.2) 

Acting on arbitrary monom ta-^ ■ ■ ■ composed from the times ta^ = tr t]""^ , where aj are 
no more restricted to be nonnegative integers, we have 
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and for the first term: 



tr |(f//,f/^ 
d 



Q2 



drjtsdrjkt 



tai ' ' ' tag 



dr]ts [fri 



tt \dvts J [\dVjJ f^jVj-Vi 



tai ' ' ' tas ■ 



(4.4) 



Thus we gave the direct proof of the "master equation" in our theory: 



2 + ^ + - aN' } Z[r^] = 0, 

where Z stands for the ratio of the determinants (G.6| 



(4.5) 



5 The genus zero solution 



Now we concentrate on the solution to Eq.(4^) in the spherical (or genus zero) hmit 
— > oo. In order to describe the solution of Eq.( [4.5| ) at large N, let us define the 
effective potential 

W[v]^ ^ log Z[r]] (5.1) 



which is normalized to be 0(1) as ^ oo. The method of solving Eq.(4^) in the large- 
limit is described in the Appendix. The result for the initial W[ri] reads 



+ br]i + c + alog (T]i + 2b + Jr]f + Ah7]i + 4c 



4Ar2 



^ + \{ri, + r]j) + c+J^ + br], + c^^ + br],+c 



(5.2) 



The variables b and c are functionals of p and a. They are determined by the nonlinear 
constraints: 

— = 0, (5.3) 



b + — y 



\ + brii + c 



and 



4A^4" 



+ c — 36^ = a. 



\ + brii + c 



(5.4) 



Due to these conditions the expression Eg. ( |5.2|) is stationary w.r.t. the differentiation 
over b and c. 



Among others, Eq. (|5.2| ) possesses a remarkable property. If we take the double deriva- 
tive of W[ri] over a, then, again due to Eqs. (|5.3|) and ( ^^) , we get 



H^M = log(62-c), 



(5.5) 



and the result coincides with that of applying the partial derivative d"^ /da^ to iy[r7]. It is 
worth to note that the answer is a pure logarithm. We shall show in a moment that the 
logarithm is immediately connected to the logarithmic scaling violation for c = 1 models. 

To deal with Eqs.( ^.3| ), (|5.4]) , let us expand the l.h.s.'s in l/rj. Using the relation 



1 



Vl — 2xy + y'^ 



n=0 



(5.6) 



where Pn{x) are the standard Legendre polynomials which are normalized by P„(l) = 1, 

(5.7) 



we rewrite Eq.( |5.3D as 




n=l 



and Eq. (p.4|) as 



where 



~ iV ^ r/f 



n> 1 . 



(5.8) 
(5.9) 



Some comments are now in order. The transformation from the variables rji to 
given by Eq.( ^.9|) is that of the Kontsevich-Miwa type for the (generalized) Kontsevich 
model. The variables become independent as — > cxd which is just our case. The 
role of to, the cosmological constant, is now played by —a. Eqs. (p^) , ( p.8|) resemble the 
corresponding equation for the Kontsevich model |§] , which is in that case nothing but the 
genus zero string equation, while we have now two variables b and c. The extra shift of t2 
in Eq.( ^.9|) is similar to that in the Kontsevich model (or in 2D topological gravity) where 
it is associated with a perturbative background. Eqs.(|5l7|), ( |5.8| ) can easily be solved in 
the following cases: 

(i) For rji ^ oo when t2 = —1, ti = t^ = . . . = 0. The solution is 6 = 0, c = a so that 

Eq. (|5.5|) gives log (—a) which recovers the known solution of the Penner model in 
genus zero. 

(ii) If all T^j's are constant: rji = s. One obtains the one-cut solution of the generalized 

Penner model with the Gaussian term added \TE, |T^. In this case Eqs.( |5.3| ) and 
lead to the relations 



— -b {-r + b+ - 

4b 2 \Ab 2 



(5.10) 
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and 



a 



-3b 



1 , s 



(5.11) 



Ab 2. 

There are two possible singular points: the first is a = resulting from 1/46 + 
b + s/2 = 0, and it is the auxiliary singularity of the Penner model. The second 
possibility is 1/46 — 6 — s/2 = 0, and it immediately gives us a = —1, without any 
reference to the value of s. It is just the true nonperturbative critical point of the 
standard Penner model. 

(iii) For a pure Gaussian case a = 0. One gets b = —\fc (remember that -Pn(l) = 1). 

One can construct perturbations of these models by finding the corresponding perturbative 
solutions of Eqs.(l^), ( p. 81 ) around the discussed ones. One should not take the double- 
scaling limit of the model ( p.l|) . Similarly to the (generalized) Kontsevich model, it 
describes already a continuum case. 

An example of such a perturbation of the Gaussian case (iii) by the Penner action can 
be obtained for a — > 0. Let b + \fc ~ a. Eq.(|5.7|) determines then a/c versus {tn}: 



, n-l 



0, 



E (2v^) 

n=l 

while Eg. ( |5.8| ) gives b + ^fc versus a: 

oo 
n=l 

(remember that P'JX) = n{n + 2)/2). Eq.(|5 



(5.12) 



n-l 



-2a 



(5.13) 

shows again a logarithmic dependence 



on a. 



Let us discuss now an interesting nonperturbative solutions to Eqs.(p.7D, (|^). It looks 
similar to the multi-critical one-cut solution of the hermitian one-matrix model with a 
polynomial potential, where for the K^^ multi-critical point one puts all tn = except 
n = 0,K. Let us take tx 7^ 0,ti = . . . = tx-i = tx+i = . . . = 0. In this case the solution 
to Eq. (|5.7|) reads 



b 



(5.14) 



where — 1 < x^^ < 1 is j*'* root of equation Pk_i{x) = (every P„(x) has exactly n 

zeros on the interval [—1,1]). Let us denote Pj 
Then for c and 6^ — c we have 



(K) 



PK{xf-\ all pf^ being nonzero. 



2a 



tKP) 



(K) 



2/K 



1] 



2a 



tKP 



[K) 



2/K 



(5.15) 



(5.16) 



One substitutes this solution into Eq. (|5.5| ) to see a logarithmic dependence of W on a. 
Therefore we conclude that our solution corresponds to a c = 1 theory. 
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In order to distinguish between different multi-critical points which are labeled by 
let us calculate the anomalous dimensions of proper conformal operators. This can be 
done by calculating their correlation functions with the puncture operator. For the partial 
derivative of W , given by Eg. ( ^.21) , over a, one gets 

— [r/] = a log(c - 6^) - c + 36^ + — ^ log [r]i + 26 + ^^i[^^^^i~+Ac J . (5.17) 

It is easy to verify that the r.h.s. is again stationary w.r.t. h and c due to Eqs.f ^^Tl) , ( ^Tsl) . 

The correlation functions can now be obtained expanding the r.h.s. of Eq. (|5.17|) in ^ 
and representing the result via t^'s. The exact answer reads 

^^Vl] = ^ E log2r/. + «log(c - h^) - f: [-^ (2v^)" . (5.18) 

We have explicitly verified this formula up to n = 4. The proof can be done by showing 
that after differentiation w.r.t. h and c and using the known rules for the derivative of 
Legendre polynomials: 

(1 - x^)PlXx) = n[P„_i(x) - xP„(x)], (5.19) 

one obtains Eqs.(^?7|), (|5.8|). 

The first term on the r.h.s. of Eq. (|5.18| ) is due to an extra factor in Eq. (|2.2|) in front 



of the integral while the remaining part looks very similar to the corresponding represen- 
tation of the derivative of the free energy of the Kontsevich model w.r.t. the cosmological 
constant in genus zero (see, e.g. |2^). Now we have two potentials h and c so that 
functions of the ratio appear which are given by Legendre polynomials. 

The correlators of proper conformal operators On with the puncture operator 1^ = 0^ 



can now be obtained by differentiating Eq. (|5.18| ). One gets 



This formula is an analog of the spherical limit of Gelfand-Dikii differential polynomials 
in the case of the Kontsevich model. For the limiting case (i) above, only correlators with 
even n = 2m are nonvanishing and one obtains 

, (2m — 1 V 



(ml 



in an agreement with an explicit calculation of Ref . |^ . 

Let us now return to our solution ( p.l5| ), ( |5.16| ). One sees that while the correlator 
of two puncture operators, given by Eq. (|0| ), depends on a logarithmically, the operator 
On is scaled as (On) ~ a"^. This dependence of n agrees with the KPZ [^| spectrum for 
c = 1. 



We think that further studies of the genus zero solutions to the model ( |2.1| ) as well as 
their extensions to higher genera deserve future investigation. 
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Appendix A 



Explicit solution at large N 



Eq.( [4.5|) can be explicitly solved as cxd while rji = 0(1). In this standard planar 

approximation all terms in the exponent ( p.4[ ) (as well as log A[a;]) are of the same order 
0{N'^). The solution (|5.2|) may be obtained along the line of the method of Brezin and 
Gross introduced for the unitary matrix model in an external field and applied to the 
hermitian one-matrix model with a cubic potential in 1^, H- 

At = cxD one defines the density of eigenvalues of t]: 



N 



p{x) 



N 



(A.l) 



i=l 



SO that W[ri\ becomes a functional of p while 



1 dW[7]\ 



N drii 



rii=x 



d 5W[p] 
dx Sp{x) 



= W{x). 



(A.2) 



Noticing that dW{x)/dx is multiplied by — the factorization property in the large-N 
limit — and can be omitted, Eg . (14.51 ) takes the form of an integral equation 



W\x) + / dyp{y) 



W {y) - W{x) 
y-x 



+ xW{x) = a for a; G [a, d] 



(A.3) 



where [a, d] is the support of p. 

It is convenient to make a shift W{x) — ^ PF(x) — |. Then Eq. (|A.3|) transforms into 



W'ix) + / dyp{y) 



W{y) - W{x) _ x^ 
y — X 4 



+ a for X G [a, d\ 



(A.4) 



where a = a + |. 

To solve Eq.(Q one reduces it to a Riemann-Hilbert problem. Let us define two 
functions 



and 



F{z) 



a Z — X 

p{x)W{x) 



dx- 



z — X 



which are analytic with cuts from a to d. It follows from Eqs.( |A.5"D ,( |A.6D that 

lm.F = W{x)\mf for xG[— oo,+oo]. 

The idea is to choose ReF to satisfy Eq. (|A.4|) , i.e. to have 



X 



YieF = W\x) + W{x)Ref - — -a 



for X G [a, d\. 

This suggests the following analytic ansatz 

F{z) = W\z) + W{z)f{z) - z^A - a, 



(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 
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where W{z) is a real analytic function restricted by 

lmW{lmW + Imf) = ioTxe[a,d] (A.IO) 

and 

lmW{Ref + 2ReW) = for x G [-00, +cx)] (A. 11) 

in order to satisfy Eqs.( |A.8D , (|A.9|) . The third restriction on W{z) is imposed by the 
asymptotic condition 

Wiz) ^^^^ z/2 + Oiz-') (A.12) 
which is a consequence of Eq. ( |A.4| ) . 

To construct W{z) that satisfies ( |A.10| )-( JA.12|) , one assumes that ImW ^ for 
X G [—/?-, where — /3+ < a to satisfy Eq.(A.lO). Let — /3_ and — /3_|_ are the roots of 

equation P'^/A + 6/5 + c = 0. Now Eq. (|A.ll|) implies 



Im I W{x)^^ + bx + cj = - bx - c for x G (A.13) 

since f{x) is real on the real axis outside [a,d]. Then Cauchy's theorem unambiguously 
determines W{x)J^ + bx + c to be 



/ 1 c b"^ (y 

+te + c= - + -61 + -- - + - 



1 /\ , ^fa) V^+fr, + c-V^ + 6. + c ^^^^^^ 

2 -^^ ./^ + 6x + c y - ^ 



where one has used Eq. ( |A.5|) and the additive polynomial has been determined to satisfy 
the asymptotic condition ( |A.12|) . 

To fix b and c we impose the requirement that W'^{z) does not have poles aX z = 
—(3± which have been assumed in obtaining Eq.( |A.ll| ) from Eq.( |A.9| ). As follows from 
Eq.( |A.14]) , the poles are eliminated providing 



1 /"^ 



p{y) 



+ by + c 



dy /^^^^ ., + c-3b' = &. (A.15) 
4 J- ^/^ + by + c 



Finally, using Eq.(A.15), W{x) can be represented as 



x2 1 r'^ p(v) x'^ + bx + c- xP^ + by + c 
W{x) = r- + bx + c+'- dy^=mJ= Vj y ^ (A.16) 



and the absence of undesirable poles is explicit in this form. 

This form of solution is also convenient to show of how the obtained W{x) satisfies 



Eq.( |A.2D . When W{x) given by ( |A.16| ) is substituted into the l.h.s. of Eq. (|A.2|) , the terms 
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which do not contain p as well as those linear and bilinear in p emerge. The key point is 
to show that the bilinear in p term splits into the product of two functionals linear in p. 
Namely, after symmetrization this term became 



dy 



p{y) 



^ + by + c 



(A.17) 



Other terms can be adjusted to eliminate the p-dependence. At this point the restrictions 
( |A.15| ) arose. 

The obtained b and c are functionals sophistically depended on p{x). To obtain W[p] 
one should integrate W{x) over x and p. The simplest way to do it is to find such 
expression for W{x) which is stationary w.r.t. the partial differentiation over b and c. 
Then we can obtain W[p] since the total derivative in p will coincide with the partial 
derivative. It appears that the expression (|A.14|) possesses the necessary property and 
after a little algebra we obtain ^^[p]: 



W[p] = f{a) + l{a-^^\ogic-b'] 



h^c -ac+^ + 3ab^ + + 
2 4 4 



+ 



rd 

/ dxp{x) 

J a 



I X^ 



- - 6 j \l — + bx + c + a\og[x + 2b + \/x^ + Abx + 4c 



dx dyp{x)p{y) log 



xy b 

— [X 

4 2^ 



(A.18) 



where the integration "constant" /(a) is not essential for applications. Similarly to the 
strong-coupling solution of Brezin and Gross 0, this W[p] is stationary w.r.t. b and c. 

Coming back to Eq.( |A.3| ), the solution to it differs from ( |A.18|) by an irrelevant to the 
critical behaviour factor —j dxx'^p{x). 
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